Berry's phase in view of quantum estimation theory, 

and 

its intrinsic relation with the complex structure 
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Abstract 

In this paper, it is pointed out that the Berry's phase is a good 
index of degree of no-commutativity of the quantum statistical model. 
Intrinsic relations between the 'complex structure' of the Hilbert space 
and Berry's phase is also discussed. 
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1 Introduction 

Berry's phase, convinced by many experiments, is naturally interpreted as 
a curvature of natural connection introduced on the line bundle over the 
space of pure states [| ||| f|l ] . 



In this paper, it is shown that the Berry's phase is a good index of non- 
commutativity of the quantum statistical model, and that Berry's phase has 
intrinsic relation with the 'complex structure' of the Hilbert space. 

The paper is organized as follows. Section [| is review of quantum es- 
timation theory and Berry's phase. In section ||, sections [5||6], relations 
between Berry's phase and quantum estimation theory are discussed. In 
section Q and sections ^-sec:timereversal, relations between Berry's phase 
and the 'complex structure' is studied. 
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2 Preliminaries 



2.1 Quantum measurement theory, the unbiased estimator 

We denote by Vi(H) the space of density operators of pure states in a 
separable Hilbert space Tt. ViiTL) is often simply denoted by V , V\. 

Let be a space of all possible outcomes of an experiment, and cr(O) be 
a a- field in f2. When the density operator of the system is p, the probability 
that the data w G lies in B € a(Q) writes 

Pi{u e B\p] = trpM (B), (1) 

by use of the map M from a (CI) to nonnegative Hermitian operator which 
satisfies 

m(4>) = o,M(n) = I, 

oo oo 

M(|J = £ M (^) ^ ^B j = cf>,i^ j), (2) 

1=1 8=1 

so that (|l]) define a probability measure (see Ref.[^], p. 53 and Ref.[10], p. 50). 
We call the map M the measurement, because there always exist an physical 



experiment corresponds to the map M which satisfies (||) P2|[17], 

The quantum estimation theory deals with identification of the density 
operator of the given physical system from the data obtained by the ex- 
periment. For simplicity, we usually assume that the density operator is a 
member of a model, or a manifold of M = {p{6)\6 £ Q C R m } C V , and 
that the finite dimensional parameter 6 is to be estimated statistically. For 
example, A4 is the set of spin states with given wave function part and un- 
known spin part. In this paper, we restrict ourselves to the pure state model 
case, where any member of the model is pure state, 

p{6) = tt(|<K0)» 

To estimate the parameter, we perform the experiment and obtain the 
data uj by which we calculate an estimator 6 by estimator 6(uj). A pair 
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(f2, er(fi), M, 9) of an estimator #(*), a measurement M, a space of data, 
a a-field cr(ri) is also called an estimator. The expectation of f(uj) with 
respect to the probability distribution (||) is denoted by Eg[f(u>)\M]. 
The estimator (fi, cx(f2), M, #) is said to be unbiased if 

E [6(u)\M\ = e (3) 

holds for all 6 € 9. If © and 

0i£ fl [0J'(u;)|M] = 6? (i,j = l,...,m), 

where <9j stands for d/d6 l , hold at a particular 0, (fl, <r(fi), M, #) is called 
locally unbiased at 0. 

As a measure of the error of the locally unbiased estimator, we use the 
weighed sum of covariance matrix, 

m 

^giMH^mu, (4) 

i=\ 

where V#[0(a;)|M] is the covariance matrix of the random variable 0{uS) which 
obeys the probability distribution (|]), and gi 9i = 1, ...,m) are non-negative 
real numbers, which corresponds to the 'cost' caused by the wrong estimation 
of the true value of the i-th component of the parameter. More generally, 

TiGV e [e(uj)\Ml (5) 

where G is a real symmetric nonnegative matrix, is also used as a measure 
of the error, and its infimum 

CR(9,G,M) = mf{TcGV e [9(u)\M] | (0, <j(0), M, 9) is locally unbiased at 9} (6) 

is called the attainable Cmnaer-Rao ( CR ) type bound, in analogy with the 
Cramer-Rao inequality in the classical estimation theory (throughout the 
paper, the term 'classical estimation theory' means the estimation theory 
of the family of probability distributions). The real symmetric nonnegative 
matrix G is called weight matrix. For CR(#, G, M), abbreviated notations 
such as CR(#, G) CR(G) are often used. 
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2.2 SLD CR inequality, the attainable CR type bound 

We have the following SLD CR inequality, which is proved for the exact 
state model by Helstrom HQ, and is proved for the pure state model by 
Fujiwara and Nagaoka ||: 

CR(0, G) > TrG ( J 5 (6>)) ~* . (7) 
Here J s (9), called SLD Fisher information matrix, is defined by 

J s (e) = \RBM9)\W))]- (8) 

where \k(0)) (i = 1, ...,m) are defined below. 

The horizontal lift \lx{9)) of a tangent vector X G T p (e)(-M) to | </>(#)), is 
an element of H which satisfies 

x P (0) = m\i x )) = \(\i x ){m\ + \m)(ix\), (9) 

and 

(i x \m) = o- (io) 

Here, X in the left hand side @ of is to be understood as a differential 
operator. |^(6>)) is defined to be a horizontal lift of di G T p ^(M) . 

The uniqueness of the horizontal lift is proved easily. As for the existence, 
given a manifold N = {\4>(0)} G 0} in 7i such that 

M = vr(AA) 

= {p(9)\p(e) = Tr(e),ee&}, (11) 

the horizontal lift explicitly writes 

|/xW) = 2X|^))-2^)X|^)). (12) 

The inequality (|^) is of special interest, because of the following theo- 
rems. 

Theorem 1 (Fujiwara and Nagaoka ffij) If the model is one dimensional, 
the equality in (|7|) establishes. 
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Theorem 2 ft] If the weight matrix is in the form of 



G = diag(0, ..,0,^,0, ...,0), (13) 
the equality in establishes. 

However, the bound is not always attainable Jl3[. 



Theorem 3 The SLD CR inequality is attainable for a strictly positive def- 
inite weight matrix, iff 

lm(l i {e)\l j {6))=0 (14) 

for any 

The model is said to be locally quasi-classical at 9 iff (|14|) holds true at 9. 
In general, 

( m 

inf \Y,9i[Ve[M]]u 



M is locally unbiased at 9 j 
> inf { gj [Vg [M] ] a \ M is locally unbiased at 9} 

i=l 

i 

holds true, where the last equality comes from theorem pi and if the model 
is not quasi-classical and the weight matrix is strictly positive, the equality 
in first inequality does not establish. The gap between the both sides of 
the inequality can be considered to be caused by the non-commutativity of 
quantum theory. 



2.3 Berry's phase 

In this section, we review the geometrical theory of Berry's phase. 

Let us denote by H the totality of state vectors, or member of 7i with 
unit length. Because the two state vectors correspond to the same state iff 
they differ with each other only in their phase factor, it is natural to consider 



5 



7i as a principal fiber bundle with the base space V\ and the structure group 
U(l). 

A horizontal lift C = {|0(t)) | < t < 1} of the curve C = {p(t) | < t < 
1} in V\ is defined to be a curve in which satisfies pit) = ir(\<fi(t))) and 

\k /d t(t)) = 2j t m))- (is) 

The connection introduced by this horizontal lift is called Pancharatnam 
connection. Then, the Berry's phase 7(C) for the curve C = {p(t) | < t < 
1} is defined by 

|0(1)) = e W>|0(l)>, (16) 

where 1 0(1)) satisfies 7r(|0(l)}) = p(l) and 

Im(0(O)|0(l))=O. (17) 

The model is said to be parallel iff Berry's phase for any curve in the model 
vanishes. 

The Berry's phase for the infinitesimal loop 

(0 1 , 9\ Oi + dP, 9 m ) <- (0\ 6 i + d6\ & + dQi, 6 m ) 

I T (18) 

e = (6\...,e\...,0i,....,e m ) -► (o 1 ,...,? + <&,...,&,...., 0™) 

is calculated up to the second order of d9 as 

-J ij dO i de j + o(d9) 2 , 

where Jjj is equal to lm(li\lj) . 
Mathematically, 

^JijdBW (19) 
is an representation of the curvature form of Pancharatnam connection. 
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3 SLD and Fubini-Study metric 



It should be noted that the SLD Fisher information matrix is deeply con- 
cerned with this fiber bundle structure. Given a curve C = {p(t)} in V±, let 
us consider the minimization 



mm ■ 



{(mm))}\<\m)) = p(t)} , (20) 



which corresponds to 'the shortest distance' between infenitesimally distant 
fibers (/?(£)) and ir~ 1 (p(t + dt)). It is pointed out that the minimum is 
achived by the horizontal lift 0, and 

For any curve C' = {\4>{t)) | < t < 1} such that C = tt(C), we can 
prove the inequality 

V t [M\t\> — ; — —. , (21) 

almost in the same way as the SLD CR inequality, where (Q, cr(fi), M, t) 
is a locally unbiased estimator of the parameter t, and dot " ' " stands for 
the differentiation with respect to t. Most strict inequality of this type is 
obtained by the minimization of the denominator of the right hand side of 
the inequality. It is already known that the minimum is given by horizontal 
lift C . Given this fact, it is easily understood that the minimum is equal to 
the SLD Fisher information matrix. 

Therefore, SLD Fisher information Jf(t) of the parameter t is propor- 
tional to 'the minimum distance' between infinitesimally distant two fibers 
ir~ 1 (p(t)) and ■K~ 1 (p(t+dt)), and the minimization corresponds to the search 
of the best possible bound. 

Remember that the inverse of the Fisher information matrix is the at- 
tainable lower bound of the covariance matrix of the unbiased estimator 
when the model is 1-dimensional. This fact implies that if SLD Fisher infor- 
mation is smaller, it is harder to distinguish p(t) and p(t + dt). We can also 
say that the closer two fibers ir~ 1 (p(t)) and 7r~ 1 ( / o(t + dt)) are, the harder 
it is to distinguish p(t) from pit + dt). 
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In the following, we define inner product (*, *}g in Tg(M) by 



(d i ,d j ) e = J s (0) 



(22) 



because this metric (SLD metric, hereafter) seems to be estimation-theoretically 
and geometrically natural. Notice the unique existence of the horizontal lift 
assures that the equality (p2l) certainly defines a metric. 



4 D-transform 

We define a linear transform D# in Tg(Ai) by 

{d il T>g{d 3 )) e = J %3 {e). (23) 

D# is called D-transform. The non-zero eigenvalues of D# are denoted by 
±i/3j(9), where (3j{0) is positive real number, and j runs from 1 to the half 
of the rank of D-transform, and (3j(0)s are sorted so that 0i > (3% > .... 
When dimM = 2, 

m = / 12(g) , (24) 

^/det J s {6) 

left hand side of which is Berry's along the curve which encloses unit area, 
where the unit of area is naturally induced by the SLD metric. 

It is worthy of remarking that the D# is a manifestation of the natural 
complex structure of the Hilbert space Ti. Actually, D-transform is obtained 
by the following procedure. 

First, multiply the imaginary unit i to \lx)- Since i\lx) is not a horizontal 
lift of any element of T p {M) generally, we project i\lx}) to spann{\li), \l m )} 
with respect to the metre Re(*|*). DJ S T p (A4) is defined to be a tangent 
vector whose horizontal lift is identical to the product of the projection. 

Defining D-transform in this way, the curevature form of Pancharatnam 
connection is defined by the equation (|23|). Therefore, Berry's phase is a 
manifestation of the complex structure of the Hilbert space. 
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multiplication of i 
\lx) G spanR{\h),..., \l m )} > i\l x ) 

T I 
taking horizontal lift project w.r.t. Re(*|*}, 

reverse image of the horizontal lift 

T i 

X e T p (M) ► Die T p (M) 

D 

5 Berry's phase and local non-commutativity 

In this section, it is pointed out that Pancharatnam curvature is good index 
of non-commutativity of two components of the parameter. As is already 
mentioned, the difference between CR(G) and TtGJ is an manifestation 
of the non-commutativity of the quantum mechanics. The author conjec- 
tures the difference increases as (3jS increase, because of the following rea- 
sons. 

First, in the 2-dimensional pure state model, we have the following the- 
orem [ 13 J . 

Theorem 4 Let Ai and Ai' be a 2-dimensional pure state model. Then, 
if/3i(0',M') > Pi(0,M) and J S (6',M') = J S (6,M), then for any weight 
matrix G, 

CR(G, 9', M') > CR(G, 6,M). (25) 

Second, by virtue of theorem ||, if Pancharatnam curvature vanishes, the 
gap between both sides of the inequality (??) vanishes. 

Third, as for general multi-dimensional pure state models, we have [13], 



CR 0,J M 



TrReJ/ m + i J 5 " 1 / 2 ^) J (9) J s -V 2 (6) 
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Notice that the transform of the parameter, 

Q i o'i = ai e\ 

and the change of the weight, 

9i = afgi, 

effects the attainable CR type bound in the exactly the same manner. There- 
fore, if every component are to be equally weighed in the sum, the weight 
matrix must be the one which 'normalize' the 'length' of each component. 
When the model is parameterized so that J s is diagonal, the weight matrix 

gives one of good normalizations, because [J s l ] u is minimum variance 
of locally unbiased estimator of the submodel of M. such that all the com- 
ponents of the parameter other than 9 1 are fixed to some constant. 

However, in the general case, the estimation theoretical meaning of 
CR(J' S ) is hard to verify. Still, this value is geometrical in the sense that 
it remains invariant under any transform of the coordinate in the model 
M. Considering the non-commutativity of the model should be invariant 
under the coordinate transform, CR(J 5 ) can be one of good measures of 
non- commut at i vity. 

6 Berry's phase and global non-commutativity 

Even if the model is locally quasi-classical at any 9 £ 0, the best locally 
unbiased estimator (f2, <r(f2), M, 9) is dependent on true value of the pa- 
rameter 9. A locally quasi-classical model is said to be quasi- classical when 
the measurement M in the best pair (Q, cr(Q), M, 9) is independent of true 
value of the parameter, for, in this case, using the globally best measurement 
M, the adaptation of the calculation of the estimate from the data, which 
is the problem of the classical estimation theory, is only left to be done to 
find the globally best pair (fi, a(Q), M, 9). 
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In the faithful model case, or the case where the model is consisted of 
the strictly positive density operator, it is pointed out the model quasi- 
classical iff Uhlmann's RPF(relative phase factor) vanishes for any curve in 
the modelKeiji:1997a. For Uhlmann's RPF is nothing but a generalization 
of Berry's phase to the non-pure model, it is of interest to study if the 
analogical fact holds true for the pure state model. 

For simplicity, we say that the manifold M in V\ is a horizontal lift of 
the model M. if 

7T (M)=M, (27) 

holds true. The horizontal lift N exists iff M is locally quasi-classical at 
any point. The model is said to be quasi-parallel iff Berry's for any curve is 
or 7T. 

Theorem 5 // the model M. is quasi-parallel, that model is quasi- classical. 



Proof First, apply Schmidt's orthonormalization to the horizontal lift M 
of A4, to obtain the orthonormal basis B = {|ej) | i = 1,2,...} such that J\f is 
a subset of the real span of B. We immerse Hilbert space H into L 2 (R, C) 

as 

I'M 61 )) = ^2 a i\ e 'i) ^ ^2ai(8)ipi(x), 

i i 

where { psiiix) \ i = 1,2,...} is an orthonormal basis in L 2 (R, C). Then, 
letting E(dx) be a projection valued measure which is obtained as a spectral 
decomposition of the position operator, the pair (R, cr(R), E, 6) is one of 
the best estimators. This assertion is easily proved by calculating the Fisher 
information matrix of the family, 

\p{x,0) = Y,(am?mx)\ 2 } (28) 
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of probability distributions. □ 

Example Let ip(x) is a wave function, or an element of L 2 (R d ,C), then, 
the model 

MMx)) = {\cp(9)) I \<f>(0)) = ip(x — 9),9 £ R d I , 

is said to be a (d- dimensional) position shifted model, and the (d- dimensional) 
momentum shifted model A4 p (ip(x)) and the (d- dimensional) position-momentum 
shifted model M x ,p(tp(x)) are defined almost in the same way. 

If ip{x) takes real value only and is symmetric about orgin (, for esample, 
the eigenstates of the harmonic oscilator), then the position shifted model 
M x (ip(x)) and the momentum shifted model M p (ip(x)) are quasi-parallel, 
and, therefore, are quasi-classical. 

Example Denote by J z the z-component of the spin operator J, by 
I to) (to = — S, — S + 1, S) the m-th eigenstate of J z , where S is the total 
spin, and define the model M. j z (\ip)) by 



MjM)=*(MjM)) 

M Jz m = urn) \ e-^\m), gr}. 

If and only if \ip) satisfies 

\{m\fP)\ = \(-m^)\ (m = -S,-S+l,...,S), 

the model Mj z {\ip)) is quasi-parallel and quasi-classical. 

Example The Riemanian Geodesic with respect to the metric tensor J s (6) 
is quasi-parallel and quasi-classical. 

The converse of the theorem is, however, not true, because the following 
counter-examples exist. 

Example We consider the model M. x which is defined by 
M x = tt{M x ) 

M x = {\<f>{9)) I \4>(d)} = cconst. x (x - ef e - {x - e)2+i9{x - e \ e e r} , 
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where g the function such that 

g{x) = 



(x > 0), 
a (x < 0). 



Then, as easily checked, N x is a horizontal lift of the model A4 X , and 
(4>{0)\4>{6')) is not real unless a = nir (n = 0,1,...). However, SLD CR 
bound is uniformly attained by the measurement obtained by the spectral 
decomposition E(dx) = \x)(x\dx of the position operator,where |xo) = S(x — 
xo). as is checked by comparing SLD Fisher information of the model J\A.x 
and the classical Fisher information of the probability distribution family 

[ P {x\9) \ P {x\e) = me)\x)\\ ogr}. 

Note that \<ft(9)} is an eigenstate of the Hamiltonian 
ft 2 d 2 K 2 ( , „ s9 1 



2m dx 2 m \ {x — 9) 2 _ 

whose potential has two wells with infinite height of wall between them. 

Example Let TL be L 2 ([0, 2tt], C), and define a one parameter model M. 
such that, 

M = ir(N) 

M= \\4>{9)) I \(/>(9)) = const, x (2-cosc;)e ia(/(LJ - 0)+9) , (0 < u, 9 < 2vr) } , 

(29) 

where a is a real number and / the function defined by 

\ u + 2vr - e (to - 6 < 0) 

Physically, ( |29"|) is an eigenstate of the Hamiltonian H such that, 

TT,n\ fi 2 ( d . \ 2 A - B cos(u; - 9) 
H(9) = la) + 



2m \du> J 2 — cos(w — 9) ' 

which characterize the dynamics of an electron confined to the one-dimensional 
ring which encircles magnetic flax $ = 2-7rac/e, where m is the mass of the 
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electron, — e the charge of the electron, c the velocity of light. A and B are 
the appropriately chosen constant. 

It is easily checked that N is a horizontal lift of the model M, and that 
the model M is not parallel unless a = nir (n = 0, 1, ...). However, consider 
the projection valued measure E w such that 

Eu(duj) = \uj) (u>\duj , 

where \loq) = S(u — u>q). Then, it is easily checked that the classical Fisher 
information of the probability distribution family 

[p{u\9) \p(u\9) = \(m\v)\ 2 , 0<u,9< 2tt] 
is equal to the SLD Fisher information of Ai. 



7 The antiunitary operator and Time reversal sym- 
metry 

7.1 The antiunitary operator 

As is pointed out in the section ||, Berry's phase seems to have some intrinsic 
relation with the 'complex structure'. In this section, we study this point 
using the antiunitary operator. 
The transformation A 

\a)=A\a), \b) = A\b) 

is said to be antiunitary iff 



<5|6> = (a\b), 
A(a\a)+0\b)) = aA\a) +pA\b), 



where z means complex conjugate of z (see Ref [18] p. 266). 

Theorem 6 The model is quasi-parallel iff the horizontal lift of the model 
is invariant by some antiunitary operator. 
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Proof Suppose that any member of the manifold N = {\4>)} in H is 
invariant by the anthmitary operator A, and let \<j)) = A\<p), \<p') = A\<p'). 
Then, we have 

<</#') = $\4>) = G R. 

Conversely, if {4>\(fi') is real for any \<j)), \(p') G M, by Schmidt's orthonormal- 
ization, we can obtain the orthonormal basis B = {\i) \i = 1,2, ...,d} such 
that M is subset of the real span of B, which means any member of N is 
invariant by the antiunitary operator which is defined by, 

K B ^2ai\i) = ^2ai\i). 

i i 

□ 

7.2 Time reversal symmetry 

As an example of the antiunitary operator, we discuss time reversal opera- 
tor (see Ref . flSfl , pp. 266-282). The time reversal operator T is an antiu- 
nitary operator in L 2 (R 3 , C) which transforms the wave function ip(x) G 
L 2 (R 3 , C) as: 

Tip(x) = ip(x) = K { \ x)} ip(x). 

The term 'time reversal' came from the fact that if ip(x,t) is a solution of 
the Schodinger equation 

then ip(x, —t) is also its solution. 

The operator T is sometimes called motion reversal operator, since it 
transforms the momentum eigenstate e* p ' X// ^ corresponding to eigenvalue p 
to the eigenstate e _ip ' X//?i corresponding to eigenvalue —p. 

Define the position shifted model by 

M x = {p(6) | p{0) = 7r(^(x - x ) ), x G R 3 }, 
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and suppose that any member of the horizontal lift A/" x of the model Ai x 
has time reversal symmetry. Then, since time reversal operator T is an- 
tiunitary, the model .M x is quasi-classical in the wider sense. The spectral 
decomposition of the position operator gives optimal measurement. 

Now, we discuss the generalization of time reversal operator. In this 
paper, the antiunitary transform 

T a : e ip yi l n — * e ia ^ e~ iv '^l n 

is also called motion reversal operator, or time reversal operator. 

If the wave function ip{x) is invariant by the time reversal operator T a , 

[ V(x - x ) V(x - x' ) dx € R (30) 

holds true for any xo, Xq, which means that the position shifted model M. x 
is quasi-parallel and, therefore, is quasi-classical. 

Conversely, if (p0|) holds true, Fourier transform of (p0|) leads to 



|*(p)| 2 = |*(-p) 

where 



(p) = — L [ V(x)e- ip ' x / ?i dx. 

V27T J 



Therefore, the wave function V( x ) is transformed to itself by the time re- 
versal operator T a such that 

T a : e ip ' X//7i — ► e i< -^ p - )+/3( ^ p - ) e~ ip '^ h 

where 



*(p) 
I*(p)I 



Theorem 7 A position shifted model is quasi-parallel if and only if there 
exists the time reversal operator which transforms the wave function ^(x) 
to itself. 
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